Fluctuation induced attraction between adhesion sites of supported membranes 
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We use scaling arguments and coarse grained Monte Carlo simulations to study the fluctuation 
mediated interactions between a pair of adhesion sites of a bilayer membrane and a supporting 
surface. We find that the potential of mean force is an infinitely long range attractive potential that 
grows logarithmically with the pair distance r: 0(r)/fcsT = clnr, where the constant c = 2 and 
c = 1 for non-stressed and stressed membranes, respectively. When, in addition to excluded volume 
repulsion, the membrane also interacts with the underlying surface through a height-dependent 
attractive potential, the potential 4>{r) is screened at large pair distances. 



Supported lipid membranes are useful and important 
model systems for studying cell membrane properties and 
membrane mediated processes 0, 0. Placing a mem- 
brane on a flat substrate allows for the application of 
several different surface sensitive techniques, including 
atomic force microscopy, x-ray and neutron diffraction, 
ellipsometry, nuclear magnetic resonance, and others 
With the aid of biochemical tools and generic engineer- 
ing, supported membranes can be functionalized with 
various membrane-associated proteins [4]. One attrac- 
tive application of supported membranes is the design 
of phantom cells exhibiting well defined adhesive proper- 
ties and receptor densities (5). Using advanced imag- 
ing techniques, detailed information can be obtained 
about the structure of the adhesion zone between the 
receptor-functionalized supported membrane and ligand- 
containing vesicles that can bind to the supported mem- 
brane [1, 0] ■ These studies provide insight into the spe- 
cific (ligand-receptor) and nonspecific interactions during 
cell adhesion [8|, |9( . Understanding these interactions is 
crucial for the development of drug delivery systems that 
depend on efficient adhesion between a liposome and the 
plasma membrane of the target cell. 

Adhesion is an immensely complex process involv- 
ing many physio-chemical and biomolecular factors [To| . 
Many aspects of this process, ranging from the cooper- 
ativity in adhesion cluster formation to the influence of 
stochastic processes such as the ligand-receptor reaction 
kinetics, have been and continue to be studied theoret- 
ically using various models (see, e.g., recent reviews in 
refs. [ll|, [ijj])- In light of this extensive theoretical ef- 
fort, it is surprising that there is still no satisfactory 
answer to one of the most fundamental problems asso- 
ciated with adhesion, namely the characterization of the 
membrane-mediated interaction between adhesion sites. 
Detailed knowledge of the strength and range of these 
interactions is essential for a better understanding of the 
role they play during the self-assembly of adhesion zones. 
The lack of theoretical studies of membrane mediated 
interactions between adhesion sites is in striking con- 
trast with the extensive literature existing on membrane- 
mediated interactions between transmembrane proteins. 
In the latter case, the origin of the interactions is the 
ability of two proteins to position themselves in a man- 
ner which minimizes the total bending elastic energy of 



the deformed membrane [llj]. In addition, the influence 
of the proteins on the membrane thermal fluctuations 
leads to membrane-mediated interactions between them 
which are analogous to Casimir forces between conduct- 
ing plates [l4| • These interactions fall off with the pro- 
tein pair separation as 1/r 4 fi1j - fil| and, therefore, at 
large distances they are considerably larger than Van der 
Waals and screened electrostatic interactions which de- 
cay much faster with r. Since adhesion sites between 
membranes or between a membrane and a surface rep- 
resent a different type of "constraint" on the shape of 
the membrane, one can expect Casimir-like interactions 
to exist between them as well. Below, we explore these 
interactions for a pair of adhesion sites between a mem- 
brane and a flat, impenetrable, surface and show that 
these interactions are of effective infinite range. 

Consider a membrane of linear size L with bending 
rigidify k and surface tension a, which may also expe- 
rience a height-dependent harmonic confining potential 
(whose second derivative with respect to the height is 7) 
due to the presence of a nearby flat surface. Let h(r) be 
the height function of the membrane, which vanishes at 
the points where the membrane is attached to the sur- 
face and takes positive vales everywhere else. The total 
energy of the membrane is, thus, is given by the effective 
Hamiltonian 



H = 



$ (h) d 2 f, (1) 



where $ represents the hard wall constraint ($ = 1 for 
h > 0, and $ = +00 for h < 0), and the integration 
is taken over the cross sectional (projected) area of the 
membranes A p ~ L 2 . Let us first consider the case where 
a = and 7 = in Eq. ([!} . In a previous publication we 
studied the behavior of a membrane with one attachment 
site to the surface Hi . We found, quite unexpectedly, 
that the attachment of the membrane to a flat surface at 
only single adhesion point does not modify the spectrum 
of thermal fluctuation of the membrane. The only effect 
of the attachment is to eliminate the membrane trans- 
lational degree of freedom by enforcing that the global 
minimum of the height function h(r) is achieved at the 
point of contact with the surface. Without the surface 
(i.e., for a freely fluctuating membrane), the manifold 
could be translated horizontally and the global minimum 



2 



could be transferred to any place within the cross sec- 
tional area r £ A p . The attachment free energy cost is, 
therefore, AF = k B Tln(A p /l 2 ) = 2k B T\n{L/l), where 
k B is Boltzmann's constant, T is the temperature, and 
I is some microscopic length scale of the order of the bi- 
layer thickness. The scaling behavior of AF with L can 
be also obtained by noting that because the single attach- 
ment point leaves the spectrum of thermal fluctuation un- 
changed, the mean height of the membrane above the sur- 
face increases as u(r) = (h(r)} ~ riy/k^T/n) with r, the 
distance from the pinning site [HI, [I?} • Helfrich showed 
that, as a result of the collisions between the membrane 
and the surface, there is an effective interaction en ergy 
per unit area: V(r) ~ (k B T) 2 / nu{r) 2 ~ (k B T)/r 2 pjj. 
By integrating this energy density over the projected area 
of the membrane, one finds that 

AF = J V(r)d 2 f= Ck B T\n ^jj . (2) 

An equation v ery similar to Eq. (f2J) has been previously 
derived in ref. [17| in the context of self-assembly of mem- 
brane junctions. In that reference, the mean field free 
energy per adhesion site was found to have a logarithmic 
dependence on the mean distance between sites. This 
result has been interpreted as a renormalization of the 
temperature downward. Our finding that C = 2, implies 
that the attachment free energy {5J exactly cancels the 
mixing entropy term of the adhesion sites and, therefore, 
the renormalizcd temperature T r = 0. The inability of 
this simple mean field calculation to predict whether the 
adhesion sites tend to aggregate (T r < 0) of segregate 
(T r > 0), emphasizes the need for a more detailed analy- 
sis of the membrane mediated interactions that also takes 
into account their many-body nature. We leave most of 
the discussion of many-body effects to a future publi- 
cation and focus here on the pair correlation function 
between adhesion sites. 

The attachment free energy ([5} is distributed within 
a volume V ~ L 2 A , where A = u(L) ~ LyJk B T/n is 
the mean height of the membrane. It, therefore, seems 
reasonable to speculate that disjoining pressure between 
the membrane and the surface scales as 



P ~ AF/V ~ y^k B ~T~K/L 3 ln(L/i). (3) 

The pressure between the membrane and the underly- 
ing surface is not uniform, however, but rather decreases 
with r because points on the manifold that are closer to 
the attachment site tend to collide more frequently with 
the surface. Defining the distance-dependent pressure, 
F(r), the mean disjoining pressure can be calculated by: 

P ~ (l/L 2 ) /j rP(r)dr, which can be reconciled with 
Eq. ([3]) by assuming that P(r) ~ ^fk B Tn/ Lr 2 . Since the 
pressure is caused by collisions between the membrane 
and the surface, one may conclude that the probability 
density that the membrane comes into contact with the 
surface at distance r from the attachment point has the 



same scaling behavior as P(r): 

n[fc(f) = o]~p(f)~i. (4) 

Let us now turn to the problem of a supported mem- 
brane with two adhesion points. Let r — fo denote 
the position of the second adhesion point within the 
cross sectional area of the membrane, while the first 
adhesion point is fixed at the origin, r = 0. The po- 
tential of mean force between the adhesion sites is de- 
fined as <j)(ro) = —k B Thi(g(ro)), where g(fo) is the 
pair distribution function expressed as a function of the 
coordinate of the second adhesion site. By definition, 
g (r ) = Z(6,r Q )/Z(6), where Z(0) and Z(0,r ) denote 
the partition functions of membranes with one (at 0) and 
two (at and fo) adhesion sites, respectively. However, 
the ratio Z(0, f )/Z(0) is also equal to n [h (r ) = 0], the 
probability density that a configuration with an adhe- 
sion site at r = makes contact with the surface at fo as 
well. We thus conclude that g(fo) = n [h (fo) = 0], and 
together with Eq. (HJ, we arrive to the following scaling 
result for the pair correlation function 

9(ro) ~ 1/rg. (5) 

From this we find that the potential of mean force be- 
tween the two adhesion sites is an infinitely long range 
attractive potential that grows logarithmically with the 
pair distance ro: 0(fo) = —k B Thi(g(ro)) = 2k B Thi(ro). 

We tested the validity of Eq. ([5]) by using constant 
surface tension (frame tension) Monte Carlo (MC) simu- 
lations of a coarse-grained implicit solvent bilayer model. 
The details of the model (which is suitable for simula- 
tions of bilayer membranes at large spatial and temporal 
scales) and simulations can be found in ref. [ill ], where 
we discuss the problem of a membrane with one attach- 
ment point to the surface. In the present work, we have 
a flat surface and 2000 coarse grained lipids that form a 
square bilayer patch of linear size L. Each lipid is rep- 
resented by a short string consisting of one head bead 
and two tail beads. The lipids reside on one side of 
the surface. Two lipids are attached to surface at their 
head beads. The location of one of these head beads 
is fixed at the origin, while the second head bead is al- 
lowed to diffuse on the flat surface. By sampling the 
position of the latter bead, the pair distribution function 
can be computed and compared with the power law dis- 
tribution, Eq. ([5]). There is, however, a problem with 
this seemingly straightforward strategy. The simulation 
time must be much longer than both (i) the typical re- 
laxation time of the longest bending mode and (ii) the 
typical diffusion time across the membrane of the mo- 
bile adhesion point. Unfortunately, both these charac- 
teristic times grow very rapidly with L, in a way which 
makes the application of the standard Metropolis MC al- 
gorithm impractical. To overcome this problem we used 
two "tricks" : The relaxation times of the thermal bend- 
ing modes were reduced by applying the recently pro- 
posed "Mode Excitation MC" (MEMC) scheme [Tg. The 
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FIG. 1: The pair correlation function, g (r), of a non-stressed 
membrane vs. the pair distance r. The slope of the dashed 
straight line is —2. 
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FIG. 2: The pair correlation function, g(r), of a stressed 
membrane vs. the pair distance r. The slopes of the solid and 
dashed straight lines are — 2 and — 1, respectively. 



MEMC scheme utilizes collective update moves that lead 
to fast excitation and relaxation of the long wavelength 
modes. The problem arising from the slow diffusion of 
the mobile adhesion point was solved by identifying the 
unpinned lipid whose head group is located closest to the 
surface and introducing a new MC move that attempts 
to place this lipid on the surface while lifting the mobile 
pinned lipid away from the surface. More precisely, if the 
height of the closest unpinned headgroup to the surface 
is ho > 0, the move attempt consists of a simultaneous 
change in the heights of all the beads of both the un- 
pinned and pinned lipids - the former are reduced by ho, 
while the latter arc increased by ho- One can easily verify 
that for this new type of MC move, detailed balance is 
satisfied by the conventional Metropolis acceptance rule: 
p(old — > new) = min(l, exp(— AE/ksT)), where AE is 
the energy change caused by the move attempt. The 
new move allows the mobile pinning site to "jump" from 
from one place on the membrane to another and enables 
efficient sampling of the pair correlation function g(fo) 
within a reasonable simulation time. 

Our results, which are shown in Fig.Q] agree very well 
with the functional form predicted in Eq. ([5]). The slope 
of the straight line on the log- log plot is equal to —2. 
Deviations from the power law behavior g(r) ~ 1/r 2 can 
be observed at small (r/L < 0.05) and large (r/L > 
0.5) distances only. At small separations, the molecular 
nature of the lipids becomes important and the radial 
pair distribution function is dominated by the depletion 
shells around the lipids. At large distances, the mobile 
adhesion site reaches the center of the square membrane 
where it becomes attracted not only by the fixed adhesion 
site but also by its periodic images. In this regime, the 
many-bony nature of the membrane-mediated potential 
must be taken into account. 

Let us now consider the case of a membrane which 
is also subjected to lateral surface tension a > (see 
Eq. ((T|)). For a stressed membrane, one can define the 
"surface tension crossover length" , £ CT ~ yj n/a, which 
marks a crossover between two regimes. On length scales 
much smaller than £ CT , the thermal fluctuations are gov- 



erned by the bending elasticity of the membrane while 
the surface tension term in Eq. ([!} is negligible. There- 
fore, for r <C one can expect Eq. ([5]) for the pair 
correlation function to hold. On length scales much 
larger than the surface tension term in Eq. ([I} be- 
comes dominant, which leads to the suppression of the 
long wavelength thermal fluctuations. Consequently, for 
r > the decay of the pair correlation function should 
be slower than predicted by Eq. ([5]). The behavior of 
g(r) in this regime can be derived using the follow- 
ing argument: Let h(r) be the height function of the 
stressed membrane. Define the function H(r) such that 

k (V 2 Lf) 2 = ct(V/i) 2 . The manifold depicted by the func- 
tion H represents a non-stressed bilayer membrane, and 
from Eq. (gj) we have that II [H (r) = 0] ~ 1/r 2 . At large 
r, the two functions can be related by the simple scal- 
ing relation kH 2 /r 4 ~ ah 2 /r 2 , i.e., h(r) — (£ CT /r)Lf(r) 
|2fj| . This scaling relationship implies that the probabil- 
ity density II [h (r) = 0] - (r/£ CT ) U[H (r) = 0]. As pre- 
viously noted, the function H[h(r) — 0] is, in fact, the 
pair distribution function, which leads us to conclude 
that in the surface-tension dominated regime, r > 
g(f) ~ (r/£er)l/r 2 ~ 1/r. The potential of mean force, 
4>(r) — — fcsTln(g(r)) = /csLln(r), simply decreases to 
half of the value of <fr(r) in non-stressed membranes. No- 
tice that this form is independent of a whose magnitude 
influences only the crossover length £ CT between the two 
scaling regimes. Fig. [2] shows our simulation results for 
a membrane with surface tension a = 3.6ksT/l 2 , where 
the microscopic length scale / is taken as the length of 
the three-bead model lipid (/ ~ 2 nm. The linear size 
of the membrane patch in our simulations L ~ 12. 51). 
The two scaling regimes with g(r) ~ 1/r 2 for small r and 
g(r) ~ 1/r for large r can be clearly seen. In a previous 
paper we measured the bending rigidity of the membrane 
and found k - 8k B T This gives \ a - 1.5/ - 0.12L, 
which is consistent with Fig. [5] where the crossover be- 
tween the scaling regimes takes place around r ~ 0.25L. 

For a non-stressed membrane experiencing a harmonic 
confining potential (7 > in Eq. (HJ)), the "harmonic con- 
finement crossover length" , £ 7 ~ (k/j) 1 / 4 , can be defined 
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FIG. 3: The pair correlation function, g(r), of a membrane 
experiencing harmonic confining potential vs. the pair dis- 
tance r. The slopes of the solid and dashed straight lines are 
—2 and 0, respectively. 

which marks the transition between two scaling regimes. 
For r -C £ 7 , the thermal fluctuations are dominated by 
the bending rigidity term in the Hamiltonian and, there- 
fore, g(r) ~ l/r 2 . For r £ T , the fluctuation spectrum 
is dominates by the harmonic confinement term. Since 
this term is a local one, the influence of the adhesion site 
is screened and pair correlation function saturates to a 
constant value, g(f) ~ r°, which means that fluctuation 
mediated force between the adhesion sites vanishes [2lj |. 
To verify these predictions, we added an energy term 
for each lipids in our model, E — (l/2)j*hf, which is 
proportional to the square of the height hi of the lipid. 
The variable 7* is related to 7 in Eq. ([1]) by 27* /a = 7, 
where a is the cross-sectional area per lipid and the fac- 
tor 2 is due to the two leaflets of the membranes. Here, 



we take 7* = 0.072, which together with the previously 
computed value a ~ 0.15Z 2 [l||, gives 7 ~ 0.096 /cgT/i 4 . 
This yields, £ 7 ~ 3^ ~ 0.25L. The results, presented in 
Fig. [3J show the two scaling regimes for small and large 
r, where the crossover between them occurs at r ~ 0.3L. 

To conclude, we have calculated the membrane- 
mediated interactions between two adhesion site of a bi- 
layer membrane and a supporting flat surface. We found 
that the potential of mean force is an infinitely long range 
attractive potential that grows logarithmically with the 
pair distance r: 0(r)/fcsT = clnr. The constant c 
takes three possible values depending on which term in 
Eq. ([1]) dominates the effective Hamiltonian of the sys- 
tem: c = 2 in the bending-rigidity dominated regime 
(which always prevails at small pair separations), c = 1 
in the surface-tension dominated regime, and c = in 
the regime where the dominant term is of the harmonic 
confinement. It is important to note that even when 
c = 2 at all pair separations (i.e., for a = and 7 = 
in Eq. Q])), the membrane- mediated attractive potential 
is not strong enough to bind the pair of adhesion sites. 
The pair correlation function in this case, g(r) ~ l/^ 2 , 
and the mean pair separation increases with the size of 
the system: (r) = J { L r 2 g(r)dr / rg[r)dr ~ L/lnZ. 
This, however, does not render the fluctuation mediated 
interactions between adhesion sites unimportant. These 
interactions may very well provide a powerful aggregation 
mechanism of receptor-ligand binding domains. In order 
to correctly analyze the aggregation behavior of an en- 
semble of adhesion sites, one must take into account the 
many-body nature of the fluctuation induced interactions 
between them. Such a study is currently underway. 
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